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MEASUREMENT  OF  POLARIZED  LIGHT  INTERACTIONS  VIA  THE  MUELLER  MATRIX 


Randall  C.  Thompson,  Jerold  R.  Bottiger,  and  Edward  S.  Fry 

Texas  A&M  University 
Department  of  Physics 
College  Station,  Texas  77843 

ABSTRACT 

A  new  instrument  for  rapid  and  accurate  measurement  of  the 
Mueller  matrix  is  described.  Distinct  measurements  of  all.  six¬ 
teen  elements  are  made  simultaneously  and  with  an  absolute  ac¬ 
curacy  of  ,1%  to  5%.  The  instrument  employs  electro-optic  modu¬ 
lators.  Results  are  presented  for  several  simple  optical  devices 


and  systems. 
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I.  INTRODUCTION 


When  light  In  some  arbitrary  polarization  state  is  incident  on 

a  medium,  the  polarization  state  of  the  scattered* transmitted,  or 

reflected  light  can  be  related  to  the  polarization  state  of  the  inci- 

1-3 

dent  light  by  a  4x4  matrix.  This  matrix,  frequently  called  the 
Mueller  matrix,  is  a  characteristic  of  the  medium.  It  depends  on 
wavelength,  and,  in  the  case  of  light  scattering,  is  a  function  of 
scattering  angle.  It  is  a  necessary  and  major  part  of  the  complete 
optical  description  of  a  medium. 

The  matrix  provides  essential  input  information  into  calcula¬ 
tions  of  the  light  field  in  planetary  atmospheres  and  in  the 
4-7 

oceans.  It  thus  has  applications  in  the  fields  of  climatology, 
illumination,  coninunications,  and  visibility.  Further  the  matri* 
itself  gives  valuable  information  about  the  biological  and  physical 

I 

properties  of  thi  scattering  particulates.  This  line  of  inquiry 


has  been  used  in 


recent  studies  of  ocean  hydrosols. 


Calculationsj 
radiation  are  in 
simple  geometries 


8-10 


of  the  matrix  for  scattering  of  electromagnetic 
|the  early  stages  of  development.  In  general,  only 
have  been  considered:  (1)  the  homogeneous  sphere 
with  arbitra/y  si  (2)  infinite  circular  cylinders;^^*^^  and 

(3)  spheroidal  particles. Recent  developments  are,  however,  showing 
considerable  promise  for  calculations  involving  arbitrarily  shaped. 


Inhomogeneous  particles 

Few  efforts  have  been  tasde  to  measure  the  entire  matrix  due  to 
the  difficulties  associated  with  using  conventiorial  techniques., 
These  involve  making  intensity  measurements  with  various  indepen¬ 
dent  combinations  of  polarizers  and  retarders  in  the  incoming  and 
outgoing  light  beams.  The  sixteen  matrix  elements  are  obtained 
from  the  frequently  small  differences  between  these  relatively  large 

measured  quantities.  Such  measurements  were  pioneered  by  Pritchard 
18  19-20 

and  Elliott  and  by  Rozenberg;  they  were  studying  the  opti¬ 
cal  properties  of  the  atmosphere.  Perhaps  the  best  example  of  the 
conventional  method  of  matrix  measurements  is  the  work  o^  Holland 
and  Gagne  ;  their  principle  objective  was  to  find  the  correlations 
between  matrix  elements  of  randomly  oriented  irregular  particles 
and  those  due  to  the  spherical  particles  having  the  "same"  size 
distribution  and  refractive  index.  Other  conventional  measurements 

directed  toward  oceanic  hydrosols  and  toward  obtaining  complete 

22 

optical  properties  of  ocean  waters  have  been  made  by  Beardsley 
and  by  Kadyshevich,  et.  al. 

Optical  polarization  modulators  can  significantly  enhance 

the  accuracy  of  polarization  measurements.  This  concept  was  first 

23 

used  in  measurements  of  the  Stokes  vector  of  a  beam  of  light. 

It  was  quite  successful  and  has  been  used  for  circular  dichroism 

24  25*26 

measurements^^  and  for  polarimetry  in  astronoiay.  "  Many  vari- 

27-30 

ants  have  appeared  in  the  recent  literature. 
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Mueller  matrix  measurements  using  polarization  modulation 

techniques  were  first  made  by  Hunt  and  Huffman^^  and  then  by 
32 

Thompson  and  Fry.  The  former  used  a  single  photoelastic  type 

modulator;  thus  they  could  in  general  only  i'«asure  combinations  of 

matrix  elements.  Nevertheless,  their  results  provided  a  ir-jor 

advance  in  the  speed  and  accuracy  with  which  matrix  trjasurements 
33 

could  be  made.  Instruments  using  a  single  modulator  have  seen 

34 

recent  use  in  biophysical  applications.  An  instrument  using  two 

polarization  modulators  has  been  recently  proposed  4nd  analyzed  by 

.  35 

Azzam. 

In  this  paper  we  will  describe  the  four  modulator  photopola* 
rimeter  developed  by  Thompson  et.  ai. 32, 36-37  provided  the 

first  direct  and  simultaneous  measurement  of  all  sixteen  elements 
of  the  Mueller  matrix.  In  light  scattering  applications  it  requires 
approximately  two  minutes  for  measurement  of  the  entire  matrix  over 
scattering  angles  from  5°  to  170®  with  an  accuracy  of  better  than 

3X.  This  compares  with  -2  hours  for  measurements  with  -lOX  accuracy 

o’. 

using  a  conventional  method.  This  ability  to  make  complete  mea¬ 
surements  very  rapidly  is  of  extreme  importance  when  one  is  studying 
samples  that  are  changing  in  time,  for  example,  ocean  water  and 
biologically  active  suspensions  In  general. 

Since  light  scattering  has  been  our  main  interest,  this  paper 
will  be  presented  from  that  point  of  view.  The  concepts,  neverthe¬ 
less,  apply  equally  well  to  transmission  and  reflection  studies. 
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II.  THEORETICAL  BACKGROUND 


To  mathematically  describe  the  po-.^Mzation  state  of  a  general 
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beam  of  light,  we  will  use  the  Stokes  vector  representation. 


S  - 


(1) 


The  linear  interaction  of  an  optical  device  with  a  beam  of  light  can 
be  described  as  a  transformation  of  an  incident  Stokes  vector  Sj 


'i 


into  an  emerging  Stokes  vector  S^, 


S,  ■  *  s, 


(2) 


Here  A  is  the  4x4  matrix  known  as  the  Mueller  (or  phase)  matrix. 

It  is  a  characteristic  of  the  optical  device  or  system.  Matrices 
representing  linear  retarders  and  polarizers  are  easily  derived.^ 
The  transformation  of  an  initial  Stokes  vector  by  a  sequence  of 
optical  devices  is  given  by  the  matrix  which  is  the  sequential  pro¬ 
duct  of  the  matrices  representing  each  device. 


A  ■  X 
-  -n 


-3^-2*  -1 


(3) 
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The  Mueller  matrix  representing  a  scattering  system  Is  given  by 
F(e),  where  e  Is  the  scattering  angle  (See  Fig.  1).  For  conven¬ 
ience  In  the  present  discussion  we  are  suppressing  the  azimuthal 
dependence;  this  is  not  a  fundamental  restriction.  For  a  system 
of  randomly  oriented  and  isotropically  distributed  particles,  F 
in  fact  has  no  azimuthal  dependence. 

The  elements  of  F  are  denoted  by  f.j(i,j»1.2,3,4).  It  is  par¬ 
ticularly  instructive  tr  work  with  the  reduced  matrix 


r(e) 


e(e) 

^11^ 


In  the  reduced  matrix  all  elements  are  restricted  by 

-1  i  i  ^  * 


(4) 


(5) 


and,  of  course,  ^il  (e)  iS' identically  plus  one.  We  will  measure 
F*(6)  directly. 

For  typical  scattering  systems  the  sixteen  elements  of  F(e)  are 

generally  not  all  independent.  Symmetries  within  the  population  of 

scattering  particles  and  within  individual  particles  serve  to  reduce 

2  1 

the  number  of  Independent  elements.  Van  de  Hulst  and  Perrin  have 
discussed  such  synmetries  in  detail.  By  measuring  the  entire  scat¬ 
tering  matrix  the  symmetries  within  it  can  be  determined  and  these 
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can  then  be  related  to  the  synmet^'Jes  In  the  scattering  system. 

Kadyshevich  et.  al.  have  found  some  symmetries  In  oceanic  hydro- 

8-10 

sols  by  using  this  approach. 
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III.  PHOTOPOLaM METER  THEORY 

In  this  section  we  win  present  the  theory  behind  the  operation 
of  our  four  modulator  photopolarimeter.  There  are,  of  course,  many 
variants  of  our  Instrument  that  eaii  be  used  to  obtain  distinct,  si¬ 
multaneous  measurement  of  all  elements  of  the  Mueller  matrix.  The 
calculatlonal  procedures  are  essentially  the  same  for  all  of  them. 
Therefore,  for  simplicity,  we  will  restrict  the  discussion  to  our 
specific  Instrument. 

A  schematic  diagram  of  uur  optical  system  Is  shown  In  Fig.  1. 

We  will  take  the  reference  plane  for  our  measurements  to  be  hori¬ 
zontal.  This  Is  Identical  to  the  scattering  plane.  MCD1,2,?, 
and  4  are  the  polarization  modulators.  They  are  electro-optic  mod¬ 
ulators  whose  optical  retardance  is  voltage  variable.  The  angle 
between  the  principle  axis  of  each  modulator  and  the  reference  plane 
Is  denoted  by  (1*1 .2,3,4).  POLl  and  P0L2  are  linear  polarizers 
with  transmission  axes  parallel  to  the  reference  plane.  The  sam¬ 
ple  whose  matrix  Is  to  be  measured  Is  denoted  by  F. 

m 

He  now  require  the  Mueller  matrix  for  each  of  the  optical  ele¬ 
ments.^  For  POLl  and  P0L2  In  the  required  orientation,  the  matrix  is 

1  1  0  Oj 

1  1  0  0 

0  0  0  0 

0  0  0  0 


(6) 


v/4  Is 


Tht  matrix  for  MODI  and  M004  with  Y^^4j  ■ 

1 
0 
0 
0 


-U4) 


For  H002  and  M003  with 

1  0 


-2(3) 


0 

d 

0 

cos 

0 

sin 

0 

1 

0 

sin  6^^^ 

0 

cos 

"2(3)  “  0 

it  is 

1 


0 

0 


0 

0 


0 

0 

cos  62(3) 

-sin  6, 

0 

0 

sin  52(3) 

cos  62 

2(3) 


(7) 


(8) 


In  the  above  matrices  (i-1 .2.3.4)  is  the  retardance  of  the  ith 
modulator.  The  matrix  to  be  measured  is  given  by 


F  • 


'll 

^12 

^13 

•'22 

^23 

'24 

'3) 

^32 

^33 

'34 

%2 

%3 

'44 

(9) 
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POLl  Insures  that  the  Initial  state  of  polarization  is  given  by 

1 

V>0  ’ 

0 


(10) 


»(her«  the  explicit  dependence  on  intensity  has  been  factored  out  as 

r.  POL  2  ensures  that  the  light  incident  on  the  PMT  has  the  Stokes 
0 

vector 


1 

1 

0 

0 


01) 


¥»here  again  the  explicit  dependence  on  the  intensity  has  been  fac¬ 
tored  out.  It  is  clear  from  Eq.  (11)  that  the  polarization  of  the 
light  incident  on  the  PMT  is  constant.  Thus  polarization  dependences 
in  the  PMT  cannot  lead  to  systematic  errors.  This  has  frequently  not 
been  the  case  in  previous  pclariza^ion  studies. 

From  Eq.  (11)  w  see  that  we  need  only  compute  the  intensity 
which  results  when  is  operated  on  by  the  optical  sys^^era.  Usitig 
the  Stokes-Mueller  calculus  we  have 


(12) 
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Using  Eqs.  6-10  and  performing  the  multiplications  indicated  in 
Eq.  (12).  we  find 


T 

^^12 

cos 

^1 

^^13 

sin 

^1 

sin 

62 

“^14 

sin 

^1 

cos 

62 

■^^21 

cos 

«4 

+f22 

cos 

«1 

cos 

«4 

■*■^23 

sin 

«1 

sin 

62 

cos 

^4 

'^24 

sin 

^1 

cos 

«2 

cos 

«4 

*^31 

sin 

*3 

sin 

«4 

*^32 

cos 

«1 

sin 

«3 

;in 

«4 

*^33 

sin 

«1 

sin 

«2 

sin 

«3 

sin 

^4 

”^34 

sin 

«1 

cos 

«2 

sin 

,'3 

sin 

'4 

"^41 

cos 

»3 

sin 

«4 

*^42 

cos 

«1 

cos 

«3 

in 

'4 

"^43 

sin 

«1 

sin 

«2 

cos 

«3 

sin 

«4 

■^44 

sin 

«1 

cos 

«2 

cos 

'3 

sin  6^) 

(13) 


All  elements  of  F  appear  Iri  this  result,  and  each  is  multiplied 
by,  a  unique  combination  of  trigonometric  functions. 

The  modulator  retardances  are  varied  sinusoidally  at  different 
frequencies 


"*1  *  *oi 


(14) 
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Here  the  amplitude  of  the  retardance  of  the  ith  modulator. 

With  this  form  for  the  retardance  the  functions  sin  5^  and  cos  5^ 
can  be  expanded  In  terms  of  Bessel  functions  of  the  retardation 
amplitudes. 

cos  S|  ■  ^“l^  *  ^*^4^ '*01^  4u^t  -  ... 

sin  4^  ■  2J^(i5q^)  cos  (i»^t  -  2J2{5q^)  cos  3u^t  +  ... 


(15) 

(16) 


These  results  are  substituted  Into  Eq.  (13)  and  all  products  are 
expanded  to  give  the  Fourier  components.  We  find  that  the  Inten¬ 
sity,  has  a  DC  component  plus  sinusoidal  terms  with  amplitude; 
proportional  to  p-oducts  of  the  Bessel  functions  and  the  matrix 
elements  f^j.  The  frequency  of  every  Fourier  component  Is  given 

k  1  lw2^inta)^  t  nw^. 
k.lfm.n*.0.1-.2.3%.... 


by 


(17) 


For  example  we  give  the  exoanslon  of  the  coefficient  of  f^^  In  equa¬ 
tion  (13), 

sin  4^  cos  42  ■  ^  "^^(^Ol^  '^0'^02^  ***1^  ”  ^  ^u^t 

"2  '12^^02^  (cos(w^  ♦  2w2)t  ^  cos  (u^  “  2ui2)t) 

♦2  J3(«o1^'^2^^02^  {cos(3w^  ♦  2w2)t  ♦  cos(3w^  «  2u.2)t) 

♦  .,..  (Id) 

Thus  the  frequencies  at  which  ^14  appears  In  the  Intensity  are  given 
by  Eq.  (17)  with  m-n»0,  h-odd  Integers,  and  l*even  Integers  and 
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xero.  The  linear  combinations  of  primary  frequencies  at  which  each 
matrix  element  appears  as  a  leading  term  In  the  Fourier  expansion 
of  the  Intensity  are  given  In  Fig.  2. 

To  Insure  the  distinct  measurement  of  each  matrix  element,  we 
will  require  that  there  be  a  unique  reference  frequency  In  the  Fourier 
spectrum  of  Eq.  (13)  corresponding  to  each  matrix  element.  Experi¬ 
mental  considerations  make  It  desirable  that  the  selected  component 
be  one  of  the  leading  terms  in  the  Fourier  expansion  and  that  the 
frequency  be  compatible  with  standard  phase  sensitive  detection 
techniques  (<100  kHz).  We  are  free  to  adjust  the  four  primary  fre¬ 
quencies,  Up  and  the  amplitude  of  the  retardances,  5Qp  sc  as  to 
satisfy  these  requirements. 

We  set  the  amplitudes  of  the  retardances  so  as  to  make 

Jq  ^^01^  ^ 

or 

-  2.404  radians.  (20), 

The  DC  term  is  thereby  removed  from  the  expansion  of  cos  6^  in  Eq. 
(15).  The  most,  important  consequence  of  this  is  that  the  DC  compo¬ 
nent  of  If  is  proportional  only  to  f^^* 

In  selecting  the  four  modulator  frequencies,  Up  we  note  that 


if  they  all  have  rational  ratios  then  the  requireinent  of  a  unique 
reference  frequency  for  each  ir.atrix  element  is  impossible  unless 
we  can  truncate  the  expansions  of  cos  and  sin  5^.  Such  a  trun¬ 
cation  is  reasonable  since  rapidly  approaches  zero  with 

increasing  n.  Furthe»*,  the  four  primary  frequencies  must  be  chosen 
so  that  the  frequencies  in  the  Fourier  expansion  of  are  not  cd 
close  as  to  require  unreasonably  long  time  constants  in  the  detec¬ 
tion  system  to  distinguish  them. 

The  set  of  primary  frequencies  we  have  chosen  for  the  present 
instrument  is 

•  41.0  kHz 
vj  ■  0.3  kHz 
vj  •  1.0  kHz 

V.  •  10.0  kHz 

'  ^  • 

The  resulting  reference  frequencies  which  we  have  selected  for  the 
sixteen  matrix  elements  are  given  in  Figure  2.  Each  of  these  re¬ 
ference  frequencies  is  indeed  unique  if  the  expansions  given  in 
Eqs.  (15)  and  (16)  are  truncated  at  Jg  and  Jg  respectively.  Extending 
each  expansion  one  more  term,  some  degeneracies  result;  the  largest 
of  these  at  our  selected  signal  frequencies  is 

^32  _  ^  «  .. 

—————  ■  . . .  3  1.9  X  10 

coeff  fj^  2  J^(«Q3)  j^(6q^) 


at  the  frequency  of  9.0  kHz.  This  is  well  within  our  exp€*'1menta1 
error  and  is  thus  acceptable.  Finally,  with  this  choice  of  primary 
frequencies,  Eqs.  (21),  all  frequencies  in  Fourier  spectrum  of 
are  separated  by  at  least  100  Hz.  Thus  they  are  easily  distinguished 
with  lock-in  detectors  and  short  time  constants. 


IV.  EXPERIMENTAL  APPARATUS 

A.  Optical  System 

The  light  source  is  a  laser  (See  Fig  We  have  used  both 

a  He-Ne  laser  (Smw.  at  633  nm)  and  a  He-Cd  laser  (20  nM  at  442  nm). 

POL  1  and  POL  2  are  Gian  prism  polarizers  which  provide  an  extinc- 
-5  ■  ■ 

tion  ratio  of  5x10  .  MODI  and  H0D2  are  the  retardance  modulators 
in  the  incident  beam.  They  are  Lasermetrics,  Model  LM-1,  trans¬ 
verse  field,  Pockels  cells.  M0D3  and  M0D4  in  the  emerging  beam 

are  longitudinal  field,  thin  plate,  Ppckels  cells  of  our  own  con- 
38 

struction.  The  crystals  in  all  four  modulators  are  potassium 
dideuterium  phosphate  (KD*P)  with  a  deuteration  level  greater  than 
99X. 


For  MODI  and  M002  approximately  750  volts  is  required  to 
achieve  the  retardation  at  633  nm.  specified  by  Eqs.  19  ai'd  20. 

For  these  modulators  the  light  must  be  collimated  and  parallel 
to  the  crystal  axisto  within  a  couple  milliradians  since  off-axis 
light  experiences  a  retardation  independent  of  the  applied  voltage. 
For  the  laser  sources,  neither  this  nor  the  relatively  small  clear 
aperture  (-2.5  tnm)  is  a  problem. 

M003  and  M0D4  generally  operate  ori  diffuse,  scattered  light. 
Therefore,  both  a  larger  clear  aperture  and  a  larger  useable'  ang¬ 
ular  field  are  required.,  Both  of  these  modulators  employ  a  single 
Z-cut  crystal  plate  12  mr  square  and  3  to  5  nw  thick.  Chromium, 


backed  gold  grids  are  vacuum  evaporated  onto  both  faces.  A  uniform 
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electric  field,  parallel  to  the  optl*:  axis,  is  produced  when  a  poten 

tial  difference  is  applied  between  the  grids.  Approximately  3000 

volts,  independent  of  crystal  thickness,  is  necessary  to  obtain  the 

required  retardation  (Eq.  20).  Since  off  axis  light  rays  undergo 

retardation  proportional  to  crystal  thickness,  the  useable  angular 

field  is  limited  by  crystal,  thickness.  However,  the  crystal  must 

be  thick  enough  to  prevent  breakdown  due  to  the  high  voltage.  This 

40 

is  particularly  a  problem  with  A.C.  voltages.  Our  experience  has 
shown  that  3  mm  and  5  tm  thicknesses  for  M003  (v^^l.O  kHz)  and  M0D4 
(10.0  kHz),  respectively,  are  good  compromises. 

The  electrc-optic  coefficient  and  the  residual  retardance  of 
KD*P  are  strongly  dependent  on  temperature.  Therefore,  to  obtain 
long  term  stability,  the  housings  of  all  four  modulators  are  held 
at  15®C  with  chilled  water  from  \  constant  temperature  (.♦0.1°C) 
circulator.  The  choice  of  15®C  is  a  compromise  between  avoidance 
of  condensation  and  reduction  of  the  required  potential  difference 
across  the  crystal. 

The  modulators  are  driven  by  a  power  amplifier  through  a  step- 
up  transformer.  For  MODI ,2,  and  3,  iron  core,  audio  transformers 
give  step  up  ratios  of  35,  80,  and  275,  respectively.  An  air  core 
transformer  giving  a  step  up  ratio  of  270  is  used  for  M0D4. 

8.  Detection  Electronics 

A  schematic  of  the  detection  system  is  shown  in  Fig.  3,  The 
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photodetector  is  an  RCA  8645  photomultiplier  tube.  A  servo  loop 
controls  the  high  voltage  to  the  photomultiplier  so  as  to  main- 

31 

tain  the  DC  component  of  the  anode  current  at  a  constant  value. 
Since  the  DC  component  is  proportional  only  to  the  effect  is 
to  divide  all  the  amplitudes  In  the  Fourier  spectrum  of  the  signal 
byf„. 

The  AC  signal  passes  through  an  AC  amplifier  and  then  to  a 
high  Q  (-200)  bandpass  amplifier  at  the  Input  of  each  of  the  fif¬ 
teen  phase  sensitive  detectors  (PSD).  The  bandpass  amplifier  is 
tuned  to  the  frequency  of  the  corresponding  matrix  element  f^j 
given  In  Fig  2. 

Reference  frequencies  for  the  PSD's  are  generated  based  on 
the  fact  that  the  product  of  two  sinusoidal  voltages  results  In 
sinusoidal  voltages  at  the  sum  and  difference  frequencies.  Ana¬ 
log  multipliers  are  used  to  form  the  products.  Bandpass  filters 
Isolate  either  the  sum  or  the  difference  frequency  and  the  re¬ 
sulting  sinusoidal  signal  Is  amplified  and  clipped  to  produce  a 
square  wave.  ■  Fig,  4  shows  a  flow;  chart  for  the  generation  of  the 
fifteen  reference  frequencies.  The  circuit  for  generation  of  a 
typical  reference  frequency  Is  shown  In  Fig.  5. 

The  phase  sensitive  detectors  (PSD)  are  based  on  analog 
multiplication  of  the  signal  from  the  photomultiplier  with  the 
appropriate  square  wave  reference  frequency.  The  Fourier  compo¬ 
nent  of  the  signal  with  the  same  frequency  as  the  reference  fre- 


quency  produces  a  DC  voltage  component  In  the  product.  This  DC 
voltage  is  proportional  to  the  matrix  element  corresponding  to  the 
reference  frequency.  The  outputs  of  all  fifteen  PSD's  (and  hence 
fifteen  matrix  elements)  are  available  simultaneously.  The  cir-  , 
cuit  for  a  typical  PSD  is  shown  in  Fig.  6. 

Only  f.|^  remains  to  be  determined.  This  can  be  done  either 
directly  in  an  independent  measurement,  or  simultaneously  by  ob¬ 
serving  the  high  voltage  to  the  photomultiplier  tube.  In  the  latter 
case,  the  servo  loop  adjusts  the  high  voltage  so  as  to  keep  the  DC 
component  of  the  signal  (and  hence  f|^)  constant.  We  empirically 
measure  this  functional  dependence  and  then  use  it  to  determine 
f-j-j  from  high  voltage  measurements.  The  high  voltage  and  the  sig¬ 
nals  from  the  fifteen  PSD’s  are  recorded  directly  on  an  analog 
recorder  or  digitized  and  stored  by  computer  for  further  analysis. 
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V.  ALIGNMENT  AND  CALIBRATION 

Proper  elignment  and  calibration  of  the  instrument  involves  the 
following  adjustments: 

a)  As  the  detector  sweeps  through  scattering  angles  it  defines 
a  plane.  The  incident  laser  beam  must  be  adjusted  so  tha^  it  lies 
in  this  plane. 

b)  The  reference  plane  for  defining  the  Stokes  vectors  must 
be  defined.  In  this  instrument  we  have  chosen  it  to  coincide  with 
the  scattering  plane.  Thus  the  transmission  axes  of  POLl  and  P0L2 
must  be  positioned  to  be  in  the  scattering  plane. 

c)  The  yaw  and  pitch  of  the  modulators  must  be  adjusted  so 
that  the  light  beam  is  parallel  to  their  optic  axes. 

d)  The  azimuth  angles  yj  of  the  modulators  must  bie  set  to  the 
values  given  in  Fig.  1. 

e)  The  retardation  amplitudes  must  be  made  to  satisfy  Eq. 

(19). 

f)  The  gains  of  the  DC  amplifiers  in  the  PSO's  must  be  set  so 

that  the  outputs  correspond  to  matrix  elements  varying  between  ±1.0. 

The  complete  alignment  and  calibration  procedure  is  described 
38 

by  Thompson.  The  adjustments  a)  and  b)  above  are  straightforward. 
We  will  briefly  describe  those  involved  in  c)  through  e). 

f 

Misalignment  of  the  light  beam  and  modulator  axes  appears  as  a 
constant  retardation.  In  a  fi^st  order  analysis  we  represent  the 
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retardance  of  each  modulator  by 

'*1  *  ^ci  ^oi  »  ^23) 

where  ‘ci  the  constant  retardation.  The  completely  general 
Mueller  matrix  representing  a  retarder  with  retardance  6  and  with 
fast  axis  at  an  arbitrary  angle  y  with  respect  to  the  reference 
plane  is  given  by  Sirmons  and  Guttmann,^ 

1  0  0  0 
0  cos^2Y+sin^2YCOsfi  h  sin4Y{l-cos6)  sin2YSin6 
0  S  sin4Y{l-cos5)  sin^2Y+cos^2YCOS5  -cos2YSin6  (24) 

jo  -sin2YSin4  coS2Y3ln6  cos4 

We  use  these  matrices  in  Eq.  (12)  to  derive  the  general  form 
of  the  detected  intensity.  In  addition  to  our  original  result  which 
assumed  perfect  alignment,  the  Fourier  expansion  of  this  result  con¬ 
tains  some  terms  at  the  four  primary  frequencies  and  some  terms  which 
mix  matrix  elements.  Specifically  the  latter  terms  are  proportional 
to  one  matrix  element  but  are  at  a  frequency  corresponding  to  a  dif¬ 
ferent  matrix  element.  These  additional  terms  are  listed  in  Fig.  7 
to  first  order  in  the  small  quantities: 

sin  ,  1  -  1,  2,  3,  4  ; 

Jq  ^^01^  .  i  -  1.  2.  3.  4  ; 

cos  (2y^)  .  i  •  1.  4  i 

sin  (2y^)  ,  i  •  2,  3  . 


(25) 
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Also,  in  Fig.  7  we  have  used  the  following  relations  which  are  cor¬ 
rect  to  first  order: 

cos  .  1  ■  U  2,  3»  4,  ; 

sin  (2y^)»1.  i  »  1.  4  ;  (zs) 

cos  (2y^)»1 ,  i  ■  2,  3  . 

Any  factor  multiplying  an  f^^  in  Fig.  7  is  to  be  made  zero  by 

the  proper  alignment  of  the  system.  The  procedure  is  to  substitute 
an  optical  device  with  known  Mueller  matrix  for  F  and  to  set  the 
detector  arm  at  e=0;  see  Fig.  1.  Specifically,  we  select  an  opti¬ 
cal  device  for  which  the  matrix  element  being  measured  is  identi¬ 
cally  zero  and  the  mixing  matrix  element  is  large.  The  appropriate 
adjustment  is  then  made  to  minimize  the  mixing  term.  The  optical 
devices  used  are  linear  polarizers,  wave  plates,  combinations  of  the 
two,  and  empty  space.  The  latter  is  the  unity  matrix. 

In  order  to  calibrate  the  full  scale  amplitudes,  optical 
devices  are  inserted  into  the  sample  space.  The  gains  of  the  ampli¬ 
fiers  in  the  PSD's  are  adjusted  to  give  2.5  volts  full  scale  if 
the  corresponding  matrix  element  is  unity.  En^ty  space  is  used 
to  calibrate  f^^,  f  33,  and  f^^.  A  rotating  crystal  polarizer  is 

used  for  f^j*  ^21’  ^13’  ^31*  ^23*  ^32’  ^  rotating  quarter-wave 

plate  is  used  for  fg^,  f^g,  f^^,  and  f^3;  the  quarter  wave  plate  re- 
tardance  is  independently  checked  by  measuring  the  previously  cali- 
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brated  Finally,  and  are  calibrated  by  using  the  quar¬ 

ter  wave  plate  and  a  rotating  polarizer  together. 

After  alignment  and  calibration,  the  accuracy  of  the  instru¬ 
ment  is  always  better  than  5%.  By  waiting  for  complete  temperature 
stabilization  and  by  carefully  iterating  the  alignment-calibration 
procedure  several  times,  the  instrument  provides  IS  measurements. 
These  accuracy  determinations  are  made  by  measuring  the  matrices 
of  known  optical  devices  and  represent  maximum  observed  errors. 
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VI.  EXPERIHENTAL  RESULTS 


Mueller  matrices  for  a  number  of  optical  devices  anti  scat¬ 
tering  systems  have  been  measured.  Some  were,  for  purposes  of  eval¬ 
uating  the  operation  of  the  instrument  while  others^^*^^  were  of 
independent  interest.  In  Fig.  8  we  shoW' the  matrix  for  a  Glan- 
Taylor  polarizer.  For  these  measurements  the  detector  is  set  at 
0°,  and  the  polarizer  is  placed  in  the  sample  space  and  rotated  at 
1.0  rpm.  The  results  are  presented  in  Fig.  8  in  matrix  form  with 
each  graph  at  the  position  of  the  corresponding  matrix  element. 

The  graphs  give  the  value  of  the  matrix  element  as  a  function  of 
the  angle  t  between  the  transmission  axis  and  the  reference  plane. 
The  normalized  matrix  for  an  ideal  linear  polarizer  is 

s1n2i|i  01 


pol 


1  cos24i 
2 

cos24)  cos  2ij» 
sin24»  HsinAii* 
0  0 


issin4*  0 
sin^2<.  0 

C  0 


(27) 


Six  of  these  matrix  elements  were,  in  fact,  used  to  calibrate  the 
instrument.  For  the  others  there  is  excellent  agreement  between 
Fig.  8  and  Eg.  (27). 

In  Fig.  9,  we  show  the  results  for  a  mounted  Mica  retarder  at 
a  wavelength  of  442  nm.  This  particular  device  had  nominally  quar¬ 
ter  wave  retardance  at  663  nm.  Data  are  obtained  in  the  same  fashion 
as  with  the  linear  polarizer.  The  abscissa  is  the  angle  y  between 
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the  retarder  fast  axis  and  the  reference  plane.  The  matrix  for 
an  ideal  linear  retarder  is  given  by  Eq.  (24).  Comparing  the  data 
with  Eq.  (24)  we  find  the  retardance  is  120°±2‘^  at  442  nm. 

Fig.  10  shows  the  Mut.'ller  matrix  for  scattering  by  a  suspen¬ 
sion  of  polystyrene  microspheres  in  water.  The  mean  diameter  of  the 
spheres  is  0.091^,  and  the  standard  deviation  is  O.OOSu.  , For -these 
data  the  abscissa  is  the  scattering  angle  9.  The  laser  wavelength 
is  633  nm.  For  this  particle  size  the  scattering  is  very  Rayleigh 
like.  Mie^^  calculations  of  the  matrix  elements  for  the  specified 
particle  distribution  are  generally  indistinguishable  from  these  data 
and  are  not  shown. 

Finally  Fig.  11  shows  the  matrix  for  a  single  polystyrene 

sphere.  The  sphere  is  suspended  in  air  in  an  inhomogeneous  electric 
42 

field.  The  laser  wavelength  is  442  nm.  Mie  calculations,  using 
l.llOu  for  the  sphere  diameter, and  1.6144  for  the  index  of  refrac¬ 
tion,  are  also  shown. 


The  agreement  is  good. 
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VII.  SUMMARY 

An  Instrument  which  provides  for  simultaneous  measurement  o^* 
all  sixteen  elements  of  Mueller  matrices  has  been  described.  Im¬ 
portant  features  are: 

(1)  A  fixed  linear  polarizer  precedes  the  photodetector.  Thus 
the  polarization  of  the  observed  light  never  changes  and  polariza¬ 
tion  sensitivities  in  the  detector  cannot  introduce  errors. 

(2)  Relatively  rapid  measurements  can  be  made  so  that  samples 
undergoing  changes  can  be  studied  (e.g.  biological  systems,  chemi¬ 
cal  reactions,  or  settling  phenomena). 

(3)  Accuracies  of  3%  are  easily  achieved  and  have  been  demon¬ 
strated  by  measuring  the  matrices  for  known  optical  devices. 

(4)  The  same  hardware  can  be  used  at  different  wavelengths 
by  simply,  adjusting  the  voltages  applied  to  the  modulators  so  as 
to  maintain  the  required  retardation  amplitudes  ag|.. 

(5)  Instrument  calibration  can  be  completely  traced  to  empty 
space  and  to  a  linear  polarizer. 

The  work  presented  here  has  been  supported  by  the  Office  of 
Naval  Research,  Code  480,  under  Grant  No,  N00014-75-C-0537. 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  of  the  instrument.  The  angle  8  is  the  scattering 
angle.  The  orientations  of  the  transmission  axes  and  fast 
axes  of  the  polarizers  and  retarders,  respectively,  are 
given  by  y. 

Fig.  2.  The  linear  combinations  of  primary  frequencies  at  which  each 
matrix  element  appears  in  first  order  in  the  output  intensity 
is  shown.  The  numerical  value. of  the  frequency  at  which  we 
observe  each  matrix  element  is  also  given. 

Fig.  3.  Schematic  of  the  electronics. 

Fig.  4.  Flow  chart  for  generation  of  the  various  reference  frequen¬ 
cies.  Analog  multipliers  are  used  to  produce  the  sum  and 
difference  frequencies. 

Fig.  5.  Circuit  for  generating  a  reference  frequency.  Standard 
operational  amplifiers  and  analog  multipliers  are  used. 
Capacitor,  resistor,  and  inductor  values  depend  on  the 
particular  reference  frequency. 

Fig.  6.  Circuit  for  phase  sensitive  detector  (PSD).  Standard 
operational  amplifiers  and  analog  multipliers  are  used. 

Fig.  7.  First  order  terms  which  arise  due  to  misalignment  of  the 
optical  system.  They  are  presented  in  the  positions  cor- 
itsponding  to  the  reference  frequencies  at  which  they  appear. 


Fig.  8.  Measured  matrix  of  a  linear  polarizer.  The  abscissa  is  the 
angle  between  the  transmission  axis  and  the  reference  plane. 

Fig.  9.  Measured  matrix  uf  a  linear  retarder  at  442  nm.  This  was 
nominally  a  quarter  wave  retarder  at  633  nm. 

Fig.  10.  Measured  matrix  of  an  aqueous  suspension. of  91.0  nm.  dia¬ 
meter  polystyrene  spheres.  The  wavelength  is  633  nm.  This 
is  a  very  Rayleigh-like  system  and  the  data  are  nearly  indis¬ 
tinguishable  from  Mie  calculations. 

Fig.  11.  Measured  matrix  for  a  single  polystyrene  sphere  suspended 

in  air  in  an  inhomogeneous  electric  field.  The  dotted  curves 
are  Mie  calculations  using  a  diamater  of  1.110  micron  and 
an  index  of  refraction  of  1.6144. 


Sch«»atlc  of  th«  Instrument.  The  .ng1«  ,  ,»  the  scetterlng 

•ngle.  The  orientations  of  the  transmission  axes  and  fast 

axes  of  the  polarixers  and  reUrders.  respectively,  are 
given  by  y. 
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Fig.  2.  The  linear  combinations  oiF  primary  frequencies  at  which  each 
matrix  element  appears  In  first  order  In  the  output  intensity 
Is  shown.  The  numerical  value  of  the  frequency  at  which  we 
observe  each  matrix  element  Is  also  given. 


Schematic  of  the  electronics 


Fig.  4.  Flow  chart  for  generation  of  the  various  reference  frequen¬ 
cies.  Analog  multipliers  are  used  to  produce  the  sum  and 
difference  frequencies* 


Circuit  for  phase  sensitive  detector  (PSD).  Standard 
operational  amplifiers  and  analog  multipliers  are  used 
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Fig.  7.  First  order  ten^s  which  arise  due  to  misalignment  cf  the 
optical  system.  They  are  presented  In  the  positions  cor¬ 
responding  to  the  reference  frequencies  at  which  they  appear. 


Measured  tiatrlx  of  a  linear  polarizer.  The  abscissa  Is  the 


near  retarder  at  442  nin.  This  was 
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tingulshable  from  Hie  calculations. 


